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ABSTRACT 

We use a combination of three-dimensional Monte Carlo radiative transfer techniques and multi- 
wavelength (BRHK, Ha) imaging data to investigate the nature of the interstellar medium (ISM) in the 
edge-on, low surface brightness (LSB) galaxy UGC 7321. Using realistic models that incorporate multiple 
scattering effects and clumping of the stars and the interstellar material, we explore the distribution and 
opacity of the interstellar material ( gas + dust ) , and its effects on the observed stellar disk luminosity 
profiles, color gradients, and rotation curve shape. We find that UGC 7321 contains a small but non- 
negligible dusty component to its ISM, yielding a B-band optical depth f e ^s ~4.0 from disk edge to 
center. A significant fraction (~ 50 ± 10%) of the interstellar material in the innermost regions of 
UGC 7321 appears to be contained in a clumpy medium, indicating that LSB galaxies can support a 
modest, multi-phase ISM structure. In spite of the clear presence of dust, we conclude that the large 
radial optical color gradients observed in UGC 7321 and other similar LSB spiral galaxies cannot be 
accounted for by dust and must result primarily from significant stellar population and/or metallicity 
gradients. We show that realistic optical depth effects will have little impact on the observed rotation 
curves of edge-on disk galaxies and cannot explain the linear, slowly rising rotation curves seen in some 
edge-on LSB spirals. Projection effects create a far larger uncertainty in recovering the true underlying 
rotation curve shape of galaxies viewed at inclinations i ;> 85°. 

Subject headings: galaxies: ISM — ISM: dust, extinction — radiative transfer — galaxies: spiral — galaxies: 
kinematics and dynamics — galaxies: individual (UGC 7321) 



1. INTRODUCTION 
1.1. Background 

Low surface brightness (LSB) disk galaxies can be 
broadly defined as rotationally-dominated galaxies with 
extrapolated face-on central surface brightness at least 
~1 magnitude fainter than the canonical Freeman (1970) 
value of fiBfi =21.65 mag arcsec -2 . Although LSB galax- 
ies span the full range of Hubble types (e.g., Schombert 
et al. 1992), here we use the term primarily in reference 
to the latest type spirals (Scd-Sdm). Numerous observa- 
tional studies have now shown that LSB galaxies are not 
"abnormal" objects, but rather are a common product of 
disk galaxy formation and evolution. More extensive sum- 
maries of the properties of LSB galaxies can be found in 
Bothun, Impey, & McGaugh (1997) and Impey & Bothun 
(1997) and references therein. 

In spite of the vast body of multiwavelength observa- 
tional data now accumulated on LSB galaxies, several as- 
pects of these systems still remain enigmatic. One im- 
portant example is the detailed nature of their interstel- 
lar medium (ISM). Although the bulk of LSB galaxies 
are known to be rich in neutral hydrogen gas (H i) , they 
are generally presumed to be dust and molecule-poor sys- 
tems. Numerous observations directly and indirectly sup- 
port this picture, including: low metal abundances (e.g., 
McGaugh 1994; Ronnback & Bergvall 1995), low statisti- 



cally inferred internal extinctions (e.g., Tully et al. 1998; 
Matthews, van Driel, & Gallagher 1998), strong similari- 
ties between optical and near-infrared (NIR) morphologies 
(Bergvall et al. 1999); undetectable CO fluxes (Schombert 
et al. 1990; Knezek 1993; de Blok & van der Hulst 1998; 
but see Matthews & Gao 2000), and the high transparency 
of their stellar disks (e.g., O'Neil et al. 1998; Matthews et 
al. 1999; Matthews, Gallagher, & van Driel 1999, here- 
after MGvD99). Theoretical and numerical models also 
suggest that the low surface densities and low metallicities 
of LSB galaxies cannot support significant molecular gas 
fractions (e.g., Mihos, Spaans, & McGaugh 1999; Gerrit- 
sen & de Blok 1999). Nonetheless, since most presently 
known LSB galaxies are actively star-forming systems, 
having blue disk colors and at least modest amounts of 
Ha emission, it seems likely that at least some molecular 
material must be present in these systems. 

Even if dust and molecular gas contents in LSB galaxies 
are typically small, extinction and reddening due to in- 
terstellar material may still have a non-negligible impact 
on the observed properties of these galaxies, and molecu- 
lar gas could still be an important constituent of the ISM 
(e.g., Spaans 1999). For example, estimates of dust red- 
dening, even if small, are needed to correctly interpret the 
stellar populations inferred from broad-band colors and 
spectroscopic measurements (e.g., Bell et al. 2000). It has 
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been found that many of the lowest-luminosity LSB galax- 
ies deviate from the standard optical Tully-Fisher relation 
in the sense that they rotate faster than predicted for their 
luminosities (e.g., Matthews, van Driel, & Gallagher 1998; 
Stil 1999; McGaugh et al. 2000), while more luminous LSB 
galaxies seem to follow the same Tully-Fisher relation as 
higher surface brightness galaxies (e.g., Zwaan et al. 1995; 
Sprayberry et al. 1995). This appears to have important 
implications for galaxy formation models (e.g., van den 
Bosch 2000) and for establishing the possible existence of 
a "baryonic" Tully-Fisher relation (e.g., Matthews, van 
Driel, & Gallagher 1998; McGaugh et al. 2000); however, 
the internal extinction corrections appropriate for the LSB 
galaxies used in these analyses remains a point of con- 
tention (e.g., Han 1992; Rhee 1996; Matthews, van Driel, 
& Gallagher 1998; Pierni 1999). Lastly, rotation curves 
of LSB galaxies offer hope to furthering our insight into 
the nature of dark matter in galaxies (e.g., de Blok & Mc- 
Gaugh 1997; McGaugh & de Blok 1998; Kravtsov et al. 
1998; Swaters, Madore, & Trewhella 2000). Accurately 
interpreting high resolution optical rotation curve data, 
which can play a crucial part in such analyses, requires 
an assessment of the role of internal extinction in these 
systems (see Section ??). In short, a more extensive and 
quantitative knowledge of the total amounts and distribu- 
tion of dust and other molecular material in LSB disks is 
overall vital to better understanding the global properties 
as well as the star-forming and evolutionary histories of 
these systems. 

To date, few attempts have been made to model the dust 
in individual LSB galaxies in a sophisticated manner. One 
difficulty is that these objects tend to be very weak far- 
infrared and sub-millimeter sources (e.g., Hoeppe et al. 
1994; Pickering & van der Hulst 1999; Bell et al. 2000), 
hence direct observations of re-emitted thermal radiation 
from their dust are lacking. Observing the dust in absorp- 
tion at optical and near-infrared wavelengths can also be 
problematic, since for LSB systems viewed at moderate or 
low inclinations, the signatures of dust can be difficult to 
infer against the patchy and diffuse background light of 
the stellar disk. The data necessary to test dust models 
are therefore sparse. For these reasons, edge-on examples 
of LSB galaxies are particularly valuable. 

Because the viewing geometries of edge-on LSB galax- 
ies provide maximal path length through the disks, such 
systems furnish us with the unique opportunity to directly 
observe the vertical and radial extent of the dust distribu- 
tion, as well as any associated dust-induced vertical and 
radial color gradients. With the aid of 3D radiative trans- 
fer models, it becomes possible to constrain the amount 
and distribution of dust and quantify its effects on the ob- 
served properties of the galaxy, including: total internal 
extinction, effects on observed color gradients, and effects 
on the disk rotation curve. Properties of the disk viewed 
at other arbitrary inclinations can also be predicted. Such 
techniques have been used to explore the effects of dust 
in more luminous, high surface brightness galaxies seen 
edge-on (e.g., Kylafis & Bahcall 1987; Kuchinski & Tern- 
drup 1996; Kuchinski et al. 1998; Trewhella, Madore, & 
Kuchinski 1999; Xilouris et al. 1997,1998,1999). Here we 
present one of the first applications of full 3-D Monte Carlo 
radiative transfer techniques to studying systems at the 
low optical depth regime of LSB galaxy disks. 



1.2. UGC 7321: A "Superthin" LSB Galaxy 

UGC 7321 is a nearby example of an Sd spiral galaxy 
with an extraordinarily thin, highly flattened and diffuse 
stellar disk with no obvious bulge component (MGvD99; 
Matthews 2000). MGvD99 showed that this "superthin" 
galaxy is actually an example of an LSB spiral galaxy 
seen near edge-on (i « 88°). MGvD99 estimated that 
after correction for projection and internal extinction, 
UGC 7321 would have a i?-band central disk surface 
brightness /jb^O) ~23.6 mag arcsec~ 2 . 

Although the observed surface brightness of UGC 7321 
is enhanced through projection, several lines of evidence 
clearly indicate this is an intrinsically LSB galaxy, in- 
dependent of internal extinction corrections. These in- 
clude: the extreme transparency of its disk (e.g., several 
background galaxies are clearly visible directly through it; 
MGvD99); the small stellar scale height and very low esti- 
mated stellar velocity dispersion (implying a very low disk 
surface density; Matthews 2000); low emission line inten- 
sity ratios (Goad & Roberts 1981); weak far-infrared and 
radio continuum emission (see below); and a rather high 
Mhi/Lb ratio (1.1 in solar units). 

UGC 7321 is an ideal LSB galaxy for detailed modelling 
since it is nearby (D « 10t| Mpc; Gallagher et al. 2000) 
and well-resolved. Some key properties of UGC 7321 are 
summarized in Table 1. Although UGC 7321 is physi- 
cally smaller, less massive, and less luminous than most 
other well-studied edge-on spirals, UGC 7321 appears to 
be a prototypical example a relatively common class of 
late-type LSB galaxies seen on edge (see also Goad & 
Roberts 1981; Karachentsev, Karachentseva, & Parnovsky 
1993; Bergvall & Ronnback 1995; Dalcanton & Schectman 
1996; Gerritsen & de Blok 1999; Matthews & van Driel 
2000). Possible face-on analogs of UGC 7321 might be 
such systems as NGC 4395 (see Sandage k Bedke 1994) 
or ESO 305-009 (see Matthews & Gallagher 1997). 

1.2.1. The Data 

Much of the primary observational data we use for our 
present analysis were described in MGvD99; these include 
near-infrared (NIR) ii-band and optical (B- & i?-band) 
imaging data, and narrow-band Ha imaging of UGC 7321. 
In addition, we supplement these data with observations 
from several additional sources. Hubble Space Telescope 
(HST) Wide Field and Planetary Camera 2 F702W and 
F814W (R- and /-band equivalent) imaging observations 
of the inner disk regions of UGC 7321 (Figure are 
an important component of our study and will be dis- 
cussed further in Gallagher et al. (2000). We also used 
the optical longslit rotation curve data of Goad & Roberts 
(1981), and NIR H and K s images of UGC 7321 obtained 
from the Two Micron All Sky Su rvey (2MASS) database 
( http: / /sirtf.jpl.nasa.gov/2mass/p . 

1.2.2. Color Gradients in UGC 7321 

One particularly intriguing finding of MGvD99 is that 
UGC 7321 exhibits quite strong radial color gradients: 
A(B — R) ~1.0 magnitude (before correcting for internal 
reddening), with the outer disk regions being considerably 
bluer than the galaxy center (Figure ||). These authors 
interpreted this as evidence of significant stellar age gra- 
dients as a function of radius in the disk, and suggested 
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that this is a galaxy that is likely to have evolved relatively 
slowly, from the inside out. Large radial color gradients 
[A(B — R) ;> 0.5] were also reported for a sample of 3 addi- 
tional superthin galaxies by Matthews, Gallagher, & van 
Driel (2000). The strength of the observed color gradients 
in UGC 7321 may also indicate that viscous evolution has 
been minimal in this disk (cf. Firmani, Hernandez, & Gal- 
lagher 1996), leading MGvD99 to postulate that galaxies 
like UGC 7321 may be some of the most pristine star- 
forming disks in the nearby universe, and hence ideal sys- 
tems in which to explore disk evolutionary processes. 

MGvD99 also reported significant vertical color gradi- 
ents in the UGC 7321 disk (as large as A(B - R) -0.45, 
with the colors reddening as a function of increasing z 
height). Such gradients are predicted to occur in galaxy 
disks as a result of dynamical heating processes (e.g., Just, 
Fuchs, & Wielen 1996), but in practice they are difficult 
to observe and interpret in most edge-on galaxies due to 
the effects of seeing, contamination from bulge light, and 
particularly due to contamination from dust. Thus rela- 
tively few empirical constraints exist for present models of 
vertical disk heating. 

A critical assumption in interpreting the physical sig- 
nificance of the large color gradients observed in the 
UGC 7321 disk, as well as the significant gradients ob- 
served in the disks of other less inclined LSB spirals (e.g., 
de Blok, van der Hulst, & Bothun 1995; Bell et al. 2000) 
is that reddening due to dust is small. For example, in a 
sample of "normal" edge-on spirals, de Grijs (1998) argued 
that to within observational errors, the observed color gra- 
dients and scale length differences in different wavebands 
could be explained solely from dust. 

Based on HST imaging observations, Matthews (1998) 
and Gallagher et al. (2000) have shown that although 
UGC 7321, as well as another nearby, edge-on LSB system 
UGC 711 lack the quintessential dust lanes seen in brighter 
edge-on spirals (cf. Howk & Savage 1999), they are not 
dust-free. Nonetheless, using a simple foreground screen 
model, MGvD99 suggested that in the case of UGC 7321, 
dust reddening could not explain more than —0.2 magni- 
tudes of the observed B — R color gradient. Using both 
foreground screen and Triplex dust models (see Disney, 
Davies, & Phillipps 1989), Bell et al. (2000) argued that 
the amount of dust reddening needed to produce the ob- 
served color gradients in their sample of less inclined LSB 
galaxies was also too large to be consistent with observa- 
tions. 

To obtain a first estimate of the role of dust in 
UGC 7321, MGvD99 employed a simple foreground screen 
model. These authors sought not to derive a realistic dust 
model for UGC 7321, but rather to test the possibility 
that a significant fraction of the observed color gradient 
in UGC 7321 could be caused by dust. To put an upper 
limit on the dust-induced color gradient, these authors as- 
sumed that the galaxy was effectively optically thin in H , 
and that the entire observed R — H radial color gradient 
could be due to dust. 

A reanalysis of the R — H color profile of UGC 7321 
(Figure ||) , actually shows a significantly larger color gra- 
dient (~ 0.9 ± 0.1 mag) than that reported by MGvD99. 
This difference was traced to an error by MGvD99 in ex- 
tracting the i?-band radial axis profile. However, because 
MGvD99 were extremely conservative in their attribution 



of the full R — H color gradient to possible dust reddening, 
as we further demonstrate here, the validity of their final 
conclusions remains unaffected. 

Although the arguments presented by MGvD99 and Bell 
et al. (2000) regarding the dust reddening in their sample 
LSB galaxies are expected to be for the most part robust, 
they nonetheless relied on simplifying assumptions, includ- 
ing the use of smooth, idealized dust geometries and the 
neglect of scattering effects. And because the color gradi- 
ents observed in UGC 7321 and some other LSB galaxies 
are often found to be large compared with those typical 
in most normal surface brightness spiral galaxies (cf. de 
Jong 1996), it is important to establish (1) whether galax- 
ies like it are in fact a genre of nearby disks with some of 
the largest intrinsic color gradients, (2) whether these are 
pathological LSBs that contain more dust than average, 
and (3) whether our knowledge of the dust in typical LSB 
galaxies is somehow incomplete. Detailed studies of indi- 
vidual galaxies rather than statistical studies are needed 
to address these questions. If the large observed color 
gradients are intrinsic (i.e., due to true stellar population 
and/or metallicity gradients) rather than a result of dust 
reddening, this is important not only for understanding the 
evolution and stellar contents of these and other LSB disks, 
but also because such gradients may imply non-negligible 
changes in the stellar mass-to- light ratio (T*) of disks as 
a function of radius. Such trends could have an impor- 
tant consequence for mass-modelling of LSB galaxies (i.e., 
for determining their dark matter halo parameters from 
rotation curves; e.g., Kent 1986) and understanding the 
underlying physics behind the Tully-Fisher relation (cf. 
Matthews, van Driel, & Gallagher 1998; Stil 1999; Mc- 
Gaugh et al. 2000). It is of interest therefore to derive 
sophisticated, realistic dust models of UGC 7321 as a first 
test case. 

We note that an advantage of studying color gradients 
in edge-on or near edge-on disks is that due to the higher 
projected surface brightnesses of the outer disk compared 
to a face-on system, we can better observe the faintest 
(and perhaps least evolved) outer disk regions. The trade- 
offs are of course the uncertain effects of dust toward the 
central disk regions, as well as the possible superposition 
of many different stellar populations along a given line of 
sight. This makes the use of full three-dimensional (3D) 
models critical for unravelling the effects of dust on the ob- 
served colors and inferred stellar population distributions. 
This is the approach we utilize in the present work. 

1.3. The Present Study 

The fundamental questions this paper seeks to address 
are: 

• What is the integrated optical depth of the dusty 
medium in UGC 7321, and what are its total mass 
and distribution? 

• What internal extinction corrections are appropriate 
for edge-on LSB galaxies? 

• How much of the observed color gradients in 
UGC 7321 and other similar galaxies can be at- 
tributed to dust? 

• What are the predicted dust-induced color gradients 
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for analogous LSB disks seen at alternate viewing 
angles? 

• What is the impact of optical depth effects on the 
shape of optically-derived rotation curves of LSB 
galaxies observed near edge-on? 

To answer these questions we use a suite of high resolu- 
tion, 3D Monte Carlo radiation transfer codes that simu- 
late scattered light images in multiple wavebands and also 
can be used to model the effects of dust on optical rota- 
tion curves. As described below, our Monte Carlo radia- 
tion transfer code has recentl y be en modified (see Wood 

to enable the efficient 
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& Reynolds 1999 and Section 
simulation of 3D, optically thin systems, thus making it 
ideal for our present study of LSB galaxies. In doing this 
we also lay the groundwork for future studies of the nature 
of dust in other LSB galaxies. Much of the general frame- 
work required for realistically modelling the effects of dust 
in galaxies has been discussed extensively by others (e.g., 
Kylafis & Bahcall 1987; Bosma et al. 1992; Xilouris et al. 
1997,1998,1999), including the application of Monte Carlo 
methods to galactic environments (e.g., de Jong 1996; Witt 
& Gordon 1996,2000; Kuchinski et al. 1998; Bianchi, Fer- 
rara, & Giovanardi 1996; Bianchi et al. 2000). We do not 
attempt to repeat the detailed findings of these works here, 
but rather apply our codes to explore specific questions 
relating to LSB galaxies. In particular, in this paper we 
focus on using our Monte Carlo models to reproduce the 
observed properties of the edge-on LSB galaxy UGC 7321 
and to predict the role dust may play on the observed 
properties of similar LSB galaxies viewed over a range of 
inclinations. 

2. THE MODELS 

The ingredients needed for our models of the dusty ISM 
of UGC 7321 are a 3D Monte Carlo radiation transfer code, 
measurements of dust and stellar distributions within the 
galaxy, and a characterization of the galactic velocity field. 
The outputs of a Mon te Carlo simulation are scattered 
light images (Sect. 3.1 & 3.2.1) and rotation curves that 
include the effects of internal extinction (Sect. ??). We 
now describe the model inputs. 

2.1. Radiation Transfer Treatment 

In the present investigation, we perform radiation trans- 
fer calculations with a 3D Monte Carlo scattering code 
that tracks photon energy packets as they are scattered 
and absorbed in a model galaxy. The code is based on 
that described by Code & Whitney (1995) and has been 
modified from the version used for the smooth density 
and emissivity models of Wood (1997) and Wood & Jones 
(1997) to enable us to study fully three dimensional sys- 
tems (Wood & Reynolds 1999). We construct the galactic 
density, emissivity, and velocity structure on a 3D linear 
Cartesian grid, and in the Monte Carlo radiation trans- 
fer, we now incorporate forced first scattering (Witt 1977) 
and a peeling off procedure (Yusef-Zadeh, Morris, & White 
1984). This enables us to very efficiently generate images 
and rotation curves for the 3D, optically thin systems we 
are investigating (see also Wood & Reynolds 1999). Unlike 
more simplistic dust treatments, our code models both ab- 
sorption and multiple scattering of photons; both are nec- 
essary for accurately determining the total amount of dust 



in a galaxy (e.g., Witt, Oliveri, & Schild 1990; Block et al. 
1994) and for predicting its observational manifestations. 

2.2. Stellar Emissivity and Dust Distributions 

In Monte Carlo simulations of the transfer of starlight 
through galaxies, the stellar sources are often represented 
by a smooth spatial distribution (e.g., Wood & Jones 1997; 
Ferrara et al. 1996, 1999; Bianchi, Ferrara, & Giovanardi 
1996) rather than by individual point sources (but see re- 
cent work by Cole, Wood, & Nordsieck 1999; Wood & 
Reynolds 1999). However, because spiral disks seldom are 
observed to have completely smooth, uniform disks at op- 
tical and NIR wavelengths, here we also consider the case 
of starlight with an additional non-uniform component su- 
perposed on a smooth background (see below). 

In a purely smooth case, for a galaxy disk that is ex- 
ponential in the radial and vertical directions, the stellar 
emissivity can be expressed as 



L(r, z) = £A,oexp 



r_ _ \z\ 
h r h z 



(1) 



where r and z are the usual cylindrical coordinates, £a,o 
is the central disk surface density at a given waveband, 
and h r and h z are the scale length and scale height, re- 
spectively, of the stars. If a bulge is present, an additional 
term will be required in the emissivity distribution. 

A smooth density distribution of interstellar material 
(i.e., gas+dust), also assumed to be exponential in both 
the radial and vertical directions, can be expressed as 



p(r,z) = p exp — 

n r ,d n z ,d 



(2) 



where po is the density (in g cm -3 ) at the galactic center, 
and h r ^ and h z ^ are the scale length and scale height, re- 
spectively, of the dust. We then define an edge-on optical 
depth in the -B-band, T ei B, measured from the disk edge 
to the disk center, as 



t e ,b = J p{r, 0)K S dr « pqKb! 



r.d 



(3) 



where the absorption coefficient kb cha ract erizes the opac- 
ity of the dust+gas mixture (see Sect. 2.3). 

Because the ISM in real galaxies is not purely smooth, 
several recent papers have examined the dust scattering of 
radiation in a two-phase medium, with emphasis on the 
penetration and escape of stellar radiation from clumpy 
environments (Boisse 1990; Witt & Gordon 1996, 2000; 
Bianchi et al. 2000). In all of these studies, clumping al- 
lowed photons to penetrate to greater depths than in a 
smooth medium, and the escape of radiation becomes en- 
hanced relative to the case where the same gas and dust 
mass was distributed smoothly. The ability of Monte Carlo 
techniques to straightforwardly incorporate these effects 
into models helps to make this technique particularly pow- 
erful, hence we also explore the role of clumping in our 
models in the present study. 

For our study of the disk of UGC 7321 we adopt the pre- 
scription for a two-phase interstellar medium from Witt & 
Gordon (1996). Here the medium consists of a smooth 
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component and a clumpy medium. This medium is char- 
acterized in terms of two parameters, namely the volume 
filling factor of dense clumps, //, and the density contrast 
between the clump and interclump medium, C. C in ef- 
fect determines the total amount of mass contained in the 
smooth versus the clumped ISM component. We empha- 
size that in the work that follows, such models generally 
contain both a clumpy and a smoothly distributed ISM 
component; however, we hereafter refer to these simply as 
"clumpy" models. 

We transform our smooth density distribution from 
Eq. 2 into a clumpy one by looping through our 3D grid 
and applying the following algorithm in each grid cell 



_/ "smooth/I// if € < //J (A\ 

Clumpy - < nsmooth /[ //( c/ - 1) + 1], otherwise. W 



where £ is a uniform random deviate in the range (0,1). 
This algorithm assures that on average the total disk mass 
is the same for a clumpy model as in a purely smooth 
model of a given optical depth, and that the ensemble av- 
erage of the clumpy distribution follows the same spatial 
profile as the smooth gas does. In this approach the small- 
est clump is a single cell in our density grid (Sect. 2.4). 

We have also applied the algorithm given by Eq. 4 to 
generate a clumpy emissivity distribution, as could result, 
for example, from clusters of young stars superposed on a 
smoother underlying disk of older stars. We discuss our 
choices for // and for C for both the stellar and gas+dust 
distributions further below. 

2.3. Dust Opacity and Scattering Properties 

We assume that the dust plus gas in our model galax- 
ies is represented by a Kim, Martin, & Hendry (1994) 
gas+dust mixture. The Kim et al. models are an exten- 
sion of the models of Mathis, Rumpl, & Nordsieck (1977) 
to include a larger distribution of grain sizes, yielding a 
higher (and more realistic) estimate of the NIR opacity. 
The wavelength dependent parameters defining this model 
are the total opacity, k\, scattering albedo, a\, and phase 
function asymmetry parameter, g\. The scattering phase 
function is approximated by the Heyney-Greenstein phase 
function (Heyney & Greenstein 1941). In Table 2 we tab- 
ulate the adopted dust parameters for our simulations. 

The total amount of dust in all of our models is charac- 
terized by the optical depth parameter r e% \ (Eq. 3), which 
is the optical depth of the model galaxy to a photon of 
wavelength A as it travels from the galaxy center through 
the disk midplane to the observer. The corresponding op- 
tical depth in the polar direction, through the galactic 
center, is 

Tp,\ = {h z .d/h rt d)T e .\. 

The Kim et al. (1994) gas+dust models that we em- 
ploy are formulated to specifically represent Galactic grain 
properties, and assumes a Galactic gas-to-dust ratio of 
140. Because the gaseous component of this ISM mixture 
contributes a negligibly small fraction of the total opacity, 
this assumption of a Galactic gas-to-dust ratio will have 
no effect on extinction or reddening effects derived from 
our models. The actual gas-to-dust ra tio ap propriate for 
UGC 7321 is discussed further in Sect. 3.6.1. 



Evidence suggests that UGC 7321 may be a rather 
metal-poor galaxy (Goad & Roberts 1981), hence it might 
be argued that Galactic grain and extinction properties 
may not be applicable, and use of an alternate extinction 
curve [e.g., that from the Small Magellanic Cloud (SMC)] 
should be considered. However, we emphasize that adop- 
tion of an SMC-type extinction curve (e.g., Bouchet et al. 
1985) would have no significant effect on the results we 
present here. The Galactic and SMC extinction curves 
are nearly identical in the optical, and begin to deviate 
significantly only in the ultraviolet, i.e. at wavelengths 
shortward of those considered in the present work (see, 
e.g., Witt & Gordon 2000). 

2.4. The Model Galaxy: Basic Structural Parameters 

To constrain the total amount of dust opacity in a galaxy 
like UGC 7321, we begin by attempting to construct mod- 
els (composed of dust, gas, and stars) that can reproduce 
the observed morphology of UGC 7321 in various wave- 
bands. Important constraints are also placed on our mod- 
els by the observed radial and vertical color gradients in 
UGC 7321, as described in detail below. 

We built our initial model galaxy based on the global 
disk parameters of UGC 7321 measured by MGvD99 and 
Matthews (2000): disk scale length h r —2.1 kpc (as mea- 
sured in R band), disk scale height h z =lA0 pc (as mea- 
sured in H band), and inclination i = 88°. For simplicity 
it was assumed that these parameters have no wavelength 
dependence. An axisymmctric, exponential distribution 
was adopted along the radial and vertical directions for 
both the dust+gas and the stars, as in Eq. 1 & 2. Using 
alternate analytic forms to describe the vertical light dis- 
tribution [e.g., a sech(z) or sech 2 (2:) function] would have 
no appreciable effect on the conclusions that follow. 

The scale height of the dust was assumed to be one half 
that of the stars {h z ,d = \h z ] e.g., Evans 1994; Xilouris 
et al. 1997,1999), and the dust scale length was taken to 
be h r> d=l-5 kpc based on the observed dust clump dis- 
tribution in the WFPC2 imaging observations shown in 
Figure || Although a number of workers have now argued 
for dust scale lengths equal to (e.g., Wainscoat, Freeman, 
& Hyland 1989) or exceeding those of the stellar light (e.g., 
Xilouris et al. 1999), such dust distributions are inconsis- 
tent with the observed distribution of dust in Figure |(see 
also MGvD99). 

The resolution of our models is set by the size of our den- 
sity grid (200 3 cells). Although the full measured radial 
extent of UGC 7321 in the i?-band is «8.1 kpc, in order to 
obtain good spatial resolution we truncated our model cal- 
culations at R m ax — ±3/i r =6.3 kpc. This yielded a spatial 
resolution of ~60 pc/pixel in our finite 200 3 cubic model 
grid and provides a good match to the projected pixel size 
of our i?-band images (~53 pc) and to the clump sizes of 
- 30 - 100 pc inferred from the WFPC2 data in Figure @ 
(see also Gallagher et al. 2000). 

Because UGC 7321 is essentially a pure disk galaxy (see 
MGvD99), no bulge component was included. Also, no 
spiral structure was added to our present models. Al- 
though we cannot know the true face-on appearance of 
UGC 7321, several arguments suggest that it is unlikely to 
exhibit well-defined spiral structure (see MGvD99). More- 
over, to first order, the effects of spiral structure are av- 
eraged over in edge-on galaxies (see Xilouris et al. 1997), 
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and Misiriotis et al. (2000) have shown via detailed models 
that derived dust and extinction parameters do not change 
appreciably if spiral structure is included. 

3. SCATTERED LIGHT IMAGES AND COLOR GRADIENTS 

As described above, we tested models for UGC 7321 us- 
ing both clumpy and smooth dust and stellar emissivity 
distributions in order to assess the amount and distribu- 
tion of dust that can best reproduce observations of the 
galaxy. In this section we examine the scattered light im- 
ages, luminosity profiles, and color gradients predicted by 
various models, and compare them to the real observa- 
tions. 

3.1. Smooth Models: Inputs and Results 

We began our modelling of UGC 7321 by generating a 
grid of models over a wide range of optical depths, and 
having purely smooth stars and dust. Because the smooth 
models are simpler (i.e., require fewer input parameters) 
than the clumpy models, these allow us to zero in on a 
suitable range of optical depths for modelling a diffuse, 
LSB galaxy like UGC 7321. Here we further examine four 
of these models, having edge-on £>-band optical depths of 
T e , B =0A, 2.0, 4.0 and 8.0 respectively. In Figure § (left), 
we show the i?-band scattered light images corresponding 
to these models. 

An examination of Figure ^ shows that we appear to 
have bracketed a suitable optical depth regime for mod- 
elling UGC 7321, but that our smooth models offer a rel- 
atively poor match to UGC 7321 morphologically. The 
models overall lack the patchy light appearance of the 
real galaxy, and in addition, the r ej s >4.0 models exhibit 
clear dust lanes in the i?-band not seen in the real data. 
Even the t e .b=2.0 model shows a brightness asymmetry 
about the midplane not seen in i?-band observations of 
UGC 7321. It is clear that in spite of the fact that an edge- 
on viewing angle averages over many of the irregularities in 
the dust and stellar distributions, a purely smooth model 
is insufficient to fully reproduce the observed properties of 
UGC 7321. Therefore in the next section consider models 
including a clumpy ISM phase, as well as a clumped stellar 
component. 

3.2. Clumpy Models: Inputs 

In this section we investigate the effects of introducing a 
clumpy, 2-phase ISM as well as a partially clumped stellar 
distribution on our Monte Carlo models. 

From our smooth models we have narrowed in on a 
range of optical depths of interest for modelling UGC 7321 
(r ej s = 0.4 — 8.0). As a next step, we must determine 
suitable values to use for the filling factor // and density 
contrast p ara meter C for both the stars and and the dust 
(see Sect. 

We note that for a clumpy medium, the optical depth, r, 
is formally somewhat ill-defined (see Witt & Gordon 2000). 
However, because our clumpy models are formulated such 
that for a given value of V they have on average the same 
disk mass as a smooth model with the same r value (see 



Sect. 2.2), one can still characterize the medium in terms 
of what is effectively a mean optical depth; we hereafter 
refer to this quantity as f . 

The clump filling factor for the ISM in real galaxies is 
very uncertain, and has been the focus of several other 



studies (e.g., Witt & Gordon 1996; Bianchi et al. 2000). As 
examples, typical values of around // = 0.15 were adopted 
by Kuchinski et al. (1998) and Witt & Gordon (2000), 
while values of 0.10-0.25 were chosen by Mihos, Spaans, & 
McGaugh (1999). However, due to the expected greater 
sparsity of molecular clouds in UGC 7321 compared with 
a higher surface brightness or Milky Way-like spiral, we 
feel the use of a somewhat lower filling factor of denser 
clumps is warranted. After trial and error, for the present 
case we find a value of ff g d = 0.06 to be suitable. 

The parameter C for the ISM is also somewhat uncer- 
tain. We therefore tested models with C=100,40,20,9, & 
2 in order span the range of values suggested by other 
workers for different environments, from Milky Way-like 
galaxies (e.g., Witt & Gordon 1996; Mihos, Spaans, & 
McGaugh 1999) to LSB galaxies (e.g., Mihos, Spaans, & 
McGaugh 1999). Our choices of C correspond to ratios of 
clumpy to smooth mass of 86:14, 72:28, 56:44, 36:64, & 
25:75, respectively. In the present case, we found the best 
qualitative matches to the data with C=20. Models with 
less than ~50% of the opacity in a clumped medium pro- 
duced image morphologies too smooth to match the data, 
while clumpier models tend to produce an overly mottled 
morphology. 

Having established estimates of the ISM parameters, we 
first ran a grid of models for various optical depths con- 
taining a clumpy ISM but purely smooth starlight. Once 
again, the resulting model images appeared too smooth 
to match the real data. Therefore these models are not 
discussed further here. Not surprisingly, we find that a 
partially clumped stellar distribution in combination with 
a multi-phase ISM will be required to realistically model 
UGC 7321. 

To model a clumped emissivity distribution, we must 
now also choose appropriate values of C» and //* for the 
stellar component of our models. In this case, we chose 
several values of C* in the range representing rough ap- 
proximations to the total amount of starlight from star- 
forming regions versus smooth underlying disk regions that 
one might expect to see in real LSB galaxies. We then ran 
grids of models with different stellar filling factors //* until 
we found good matches to the observed image morpholo- 
gies of UGC 7321 in the B, R, and H bands. We viewed 
each of our models from both an i = 88° and a face-on ori- 
entation to check that our model looked qualitatively like 
a real late-type LSB galaxy disk from multiple viewing an- 
gles. Through this scheme we arrived at preferred values 
of //*=0.60 and C*=2. This results in an apportionment 
of ^75% of the light into higher density component and 
^25% into a smooth background component. However, we 
emphasize that none of our key conclusions regarding the 
total amount of dust in UGC 7321 or its corresponding ef- 
fects on the observed properties of the galaxy are strongly 
sensitive to the choice of these values. 

Lastly, we note that one might expect some small wave- 
length dependence on C* and //* to account for the fact 
that generally the starlight in galaxies becomes intrinsi- 
cally smoother as one moves from the optical to the NIR 
regime. However, as shown by Bergvall et al. (1999), LSB 
galaxies tend to exhibit relatively minor morphological 
differences between the optical and NIR compared with 
brighter spirals. Therefore we have settled on one value of 
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//* and of C* that appear to give satisfactory matches to 
the data over the range of wavelengths considered here. 

3.2.1. Scattered Light Images from the Clumpy Models 

i?-band scattered light images for our grid of clumpy 
models with f e ,B=0.4, 2.0, 4.0, and 8.0 are shown in Fig- 
ure H (right), alongside the smooth models from Sect. [0] 
for comparison. It is immediately clear that unlike the 
smooth models, our clumpy models do not exhibit dust 
lanes within this range of optical depths. Rather one 
sees the signatures of clumpy, irregular absorption over 
the central few kpc of the galaxy, consistent with the real 
UGC 7321. All models also show a much less pronounced 
asymmetry about the midplane compared with the smooth 
models of the corresponding optical depth. In addition to 
the models shown in Figure |], we also ran clumpy models 
with f e) B=12 and 24, but these models showed clear dust 
lanes and could be immediately ruled out as having too 
much opacity. 

Of the models shown in Figure |[ the f e s=4.0 case ap- 
pears to provide the best fit to multiwavelength observa- 
tions of UGC 7321 on a purely morphological basis. In 
the f ej B=8.0 model, we see too few bright emitting re- 
gions more than one pixel across, and somewhat too much 
central opacity. In contrast, the f Bi B = 0.2 and 2.0 models 
both show stronger central concentrations of starlight at 
i?-band than in the real galaxy, and less patchiness within 
the inner few kpc. The r ei s=4.0 has a discernible but 
weak central concentration of starlight, and little evidence 
for absorption outside roughly ±3 kpc, both in excellent 
agreement with UGC 7321. The match of the f ej B=4.0 
model is also quite good at iJ-band, where it nicely re- 
produces the more well-defined brightness center of the 
galaxy and overall smoother appearance of the central re- 
gions of the galaxy as seen at these wavelengths. In Fig- 
ure H we show more detailed close-ups of the f ei B=4.0 R 
and H scattered light images compared with view of the 
real galaxy at the corresponding wavelengths. 

3.3. Color Gradients Predicted from Clumpy Models 

We now explore the color gradients predicted by our 
models and compare them to those seen in the real data. 
In Figure ^| we plot the B — R and R—H radial color gradi- 
ents along the maj or a xis extracted from both our smooth 
models (from Sect. 3.1) and our clumpy models, for the op- 
tical depths r ej B=0.4, 2.0, 4.0, and 8.0. As expected, the 
color profiles extracted from the clumpy models appear 
somewhat "noisier" than those from the smooth models. 
In fact, the intrinsic rms noise inherent in our models in 
only ^0.005 magnitudes, hence the visible fluctuations in 
the radial color profiles in the lower portion of Figure || 
are real. 

From Figure ^ we see that at a given optical depth, 
the total dust-induced color gradient (from disk edge to 
center) is virtually identical to that in the smooth mod- 
els. Likewise, the mean amount of reddening at a given 
galactocentric radius is nearly identical in the two sets of 
models; along a given sight line, differences appear to be 
<0.015 magnitude. In general, the reddening effect of a 
clumpy medium is expected to be smaller than that pro- 
duced by a diffuse model with the same amount of dust, 
particularly if the filling factor of the clouds is small (e.g., 
de Jong 1996; Bianchi et al. 2000). However, our new 



models show that in the low optical depth regime, even 
for the low clump filling factors used here, this effect is 
nearly imperceptibly small, and that both the clumpy and 
the smooth models produce roughly the same amount of 
mean reddening, albeit with much larger statistical varia- 
tions between different lines of sight. We find this holds 
even for a pure clumped medium with no smooth compo- 
nent. 

A further examination of Figure ^ reveals that the ob- 
served A(R — H) ^0.9 of Figure || cannot be reproduced 
solely from dust effects for models with t &: b < 8.0. Mean- 
while, the B — R color profiles for both the smooth and the 
clumpy models over this optical depth range show an addi- 
tional interesting behavior: for r 6) B >4.0 the dust-induced 
B — R color actually saturates at A(B — R) ~0.31. In other 
words, adding additional dust cannot further increase the 
total dust-induced B — R color gradient to >0.31 magni- 
tudes from the disk edge to the disk center. As seen from 
Figure ||, A(B — R) ^0.3 is significantly less than the ac- 
tual observed B — R color gradient in UGC 7321 over this 
radial interval, hence our models indicate that no amount 
of dust can fully account for the observed B — R color 
gradient in UGC 7321. 

This saturation effect we see in the dust-induced B — R 
color gradients arises due to the wavelength-dependent ab- 
sorption in realistic systems where emission and absorp- 
tion are mixed together (in contrast, for example, to simple 
foreground screen models). A simple illustrative example 
is a constant density absorbing medium (no scattering) 
with uniform emissivity, i.e., a slab model where the stars 
and dust are spatially coincident. The emergent intensity 
is I\ oc [1 — e~ Tx ]/r\. From Table 2, we see that (approx- 
imately) Tb • Tr : Th = 8:4:1, giving 



Jb [l-e~ r *] Ir ^ [1 - e-^/ 2 ] 
7b X 2[1 - e-W 2 ] ' I H ° C 4[1 - e - T «/ 8 ] 



(5) 



Plotting (B — R) and (R — H) against tb for this simple 
example (Figure 0), we see that (B — R) saturates be- 
fore (i? — H) due to the smaller opacity difference in going 
from B to R, than in going from R to H. The exact optical 
depth at which the saturation occurs is slightly different 
for this simple model than the more realistic models we 
have computed. However it still demonstrates that even if 
we have significantly underestimated the amount of opac- 
ity in UGC 7321's disk, our conclusion remains unchanged 
that dust alone still cannot account for the large color gra- 
dients observed in UGC 7321. 

Above we have argued that models for UGC 7321 with 
T e ,B >8.0 can be ruled out on a purely morphological basis. 
To further demonstrate that we have not grossly underes- 
timated the amount of internal extinction in UGC 7321, 
we have also examined the H — K color gradients. In Fig- 
ure ^ we show the H — K color profile along the major axis 
for our family of clumpy models, and in Figure |^ we show 
a radial H — K s color plot for UGC 7321 derived from 
2MASS data. The field-of-view of the available 2MASS 
data is limited to ~ 200", and due to the large angular 
size of UGC 7321 compared with the scan size of the sur- 
vey, the data have very poor signal-to-noise outside the 
central arcminute or so of the galaxy (see Jarrett et al. 
2000 for an description of 2MASS extended source data). 
Nonetheless, the 2MASS data are consistent with a H — K 
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radial color gradient of no more than ~0.1 magnitude in 
UGC 7321. Comparing this to our models, we see that 
f e ,B=S.O models produce a A(H — K) in excess of that 
expected on the basis of the 2MASS results. 

As a final test on our estimates of r, we lastly examine 
the vertical B — R color gradients (i.e., those parallel to 
the minor axis) for the same four optical depths as above 
(Fi gur e [nj) and compare these with the real data (Fig- 
ure O). 

In the real UGC 7321 data, near the minor axis, the 
vertical B — R color profile is seen to be redder by —0.31 
magnitudes at z=0.15 kpc compared to the highest ob- 
servable z-heights (Figure [ll]). However, we see that the 
T e ,B=8-0 model predicts a reddening of several tenths of a 
magnitude higher than this along the minor axis. Mean- 
while, since reddening toward the midplanc along the z 
direction is in general due to dust rather than stellar pop- 
ulation effects, it appears f e ,s >2.0 is required in order to 
explain the amount of reddening that is observed. Hence 
we again conclude that the f ei s=4.0 model produces the 
best overall match to the data. 

Near a galactocentric radius of V.O (~3 kpc), the f e< B — 
4.0 model predicts a dust-induced vertical B — R color gra- 
dient of —0.15 magnitudes. In the real UGC 7321 data, 
MGvD99 actually observed a bluing of A(B - R) -0.3 
toward the galaxy midplane at this radius. Such a bluing 
toward low z-heights is predicted to occur as a consequence 
of new young stars continually being born near the galaxy 
midplane, while older stars are scattered to larger z heights 
via dynamical heating processes (e.g., Just, Fuchs, & Wie- 
len 1996). The small amount of dust reddening helps to 
explain why UGC 7321 is one of the few galaxies studied 
so far in which we can readily observe the predicted stellar 
population changes with z-height. This demonstrates that 
edge-on LSB galaxies like UGC 7321 can provide excel- 
lent laboratories for further constraining dynamical heat- 
ing processes. 

In summary, we conclude that a clumpy ISM model with 
T e ,B ~ 4-0 provides the best characterization of the observed 
properties of UGC 7321, and that dust alone cannot ac- 
count for the large radial color gradients observed in this 
galaxy. This suggests that the existence of significant stel- 
lar population gradients (and possibly metallicity gradi- 
ents) in the disk of UGC 7321, and the likelihood that 
galaxies like UGC 7321 maintain a well-preserved record 
of how galaxy disks have been built up over time. 

3.4. Radial Intensity Profiles 

In Figure [l2] we show B, R, and H radial brightness 
profiles extracted along the major axes for our grid of 
clumpy models. These profiles permit estimates of the 
expected extinction along different lines of sight and in 
different wavebands for a galaxy like UGC 7321. These 
profiles confirm that for a f ej s=4.0 model, extinction ef- 
fects are nearly negligible in iJ-band, even at small galac- 
tocentric radii (see also MGvD99; Matthews 2000). They 
also show that in B-band, extinction toward the the center 
of UGC 7321 should be roughly 0.65±0.2 magnitudes in 
B-band and 0.43±0.2 magnitudes in R, in good agreement 
with the estimates made by MGvD99. 

3.5. UGC 7321 Analogs at Other Inclinations 



We have seen from our models in the previous sections 
that UGC 7321 is dust poor, but not dust free, and that 
small but non-negligible corrections must be made for the 
total internal extinction and reddening caused by dust 
when interpreting the colors, total magnitudes, and in- 
trinsic surface brightnesses of LSB galaxies viewed at high 
inclination. But what about at lower inclinations? 

An advantage of our 3D Monte Carlo simulations is the 
ability to easily examine a given model at any arbitrary 
viewing angle. We can therefore use the models described 
above to predict how dust will affect the observed prop- 
erties of LSB disk galaxies seen at inclinations higher or 
lower than 88° . Other authors have previously used Monte 
Carlo simulations to explore the effects of dust on the col- 
ors and global extinctions of less inclined galaxies (e.g., de 
Jong 1996), but our high resolution, high signal-to- noise 
simulations allow us to probe the detailed effects of dust 
on the observed light and color distributions in the very 
low optical depth regime appropriate for LSB galaxies and 
also to investigate the effects of a clumpy ISM. All of the 
models we discuss in this section have stellar and dust 
scale parameters identical to those used in the UGC 7321 
models of the previous sections (see also Table 2). 

In Figure [l3| & [l4| we show the predicted R—H and B—R 
major axis color profiles for our r Ei B=0. 4-8.0 smooth and 
clumpy models, respectively, at inclinations i = 0°, 45°, 
and 90°. For comparison, we also replot our i = 88° mod- 
els from Figure [| 

At i = 90°, the observed trends are much as expected. 
For the various optical depths, we see small increases of or- 
der 0.1 magnitude in the amplitude of the R — H color gra- 
dients compared with the i = 88° models, and difference of 
— 0.05 magnitude in A(B — R) compared with the i = 88° 
case. The effects at i = and i = 45° are somewhat more 
interesting, although primarily from a theoretical perspec- 
tive. In cases of galaxies viewed at high inclinations, light 
scattered by dust has only a small effect on the emergent 
light profiles, particularly in a clumpy medium. Instead 
the observed color changes are dominated by dust absorp- 
tion effects. In contrast, for i < 45° and the small opti- 
cal depths we consider here, we can see directly from our 
models how the dominance of scattering versus absorption 
changes across the face of the galaxy. In these cases, dust 
scattering actually causes a bluing of the light towards the 
galaxy center. For the higher optical depth (t &: b >2.0) 
models, this trend is interrupted by a rather abrupt red- 
dening at the smallest galactocentric radii. These trends 
are easier to see in the smooth models (Figure |l3|). How- 
ever, note that at both i = and i = 45°, the maximum 
amplitude of these effects are only a few hundredths of a 
magnitude and therefore are predicted to essentially un- 
observable in real LSB galaxies. We conclude that at low 
and moderate optical depths, dust will have no appreciable 
effect on the observed radial color gradients of moderately 
inclined galaxy disks, consistent with the findings of de 
Jong (1996). Likewise, internal extinction corrections at 
optical and NIR wavelengths should be almost negligible 
for the bulk of LSB galaxies viewed at low or moderate 
inclinations. 

3.6. Discussion 
3.6.1. Gas-to-Dust Ratios and Total Dust Content 
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In Sections 3.2.1 & 3.3, we concluded that among the 
models we have considered, a clumpy, multi-phase ISM 
model with f ej s w 4.0 provides the best overall match 
to the observed properties of the edge-on LSB galaxy 
UGC 7321. We now consider some additional implications 
of this model. 



As described in Sect. 2.2, we have parameterized our 
models in terms of the opacity of a gas+dust mixture. For 
a given optical depth, the total mass of material in our 
model ISM can then be computed as 



M ISM = / pdV 



2np Q 



exp 



h z ,d 



(6) 
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From this, one can obtain the mass in the form of dust 
by dividing Eq. 6 by the assumed gas-to-dust ratio of the 
absorbing medium (in our case, p g / p d —lA0; see Sect. 2.3). 
In terms of the edge-on optical depth t e ,B , the dust mass 
M d can thus be expressed in solar units as 



M d = 2.2 x 10 7 • (pg/pd) T~e,B M G 



(7) 



For our f e s=4.0 clumpy models, this yields M d ~ 8.8 x 
10 5 M.q, with 56% of this in a clumpy component, and 
the remainder in a smoothly distributed medium. 

IRAS far-infrared flux measurements are also available 
for UGC 7321 (e.g., Sage 1993), allowing an independent 
estimate of the dust content of UGC 7321. From the IRAS 
100/im flux, one can make an estimate of the mass of the 
far-infrared-emitting dust in a galaxy as: 
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where Sioo is the IRAS 100/im flux density in Jy (1.15 Jy 
for UGC 7321; Sage 1993), D is the distance in Mpc, T d 
is the dust temperature, and c g is a constant that de- 
pends on the grain opacity, and can vary by roughly a 
factor of two, depending on the composition of the grains 
(e.g., Devereux & Young 1990). Sage (1993) estimated 
7^=33.5 K for UGC 7321, and we adopt the value of 
c g =5 suggested by Thronson & Telesco (1986). This yields 
M dtWOlim ~ 4.2 x 10 4 M Q for UGC 7321, or only -5% of 
the dust mass estimated from our models. However, De- 
vereux & Young (1990) found from an examination of a 
large sample of JiMS'-detected nearby galaxies, that typi- 
cally only ^ 10-20% of the total dust mass in normal disk 
galaxies is warm enough to radiate in the IRAS bands (see 
also Greenberg & Li 1991; Whittet 1992). Lisenfeld & Fer- 
rara (1998) also found for a sample of dwarf galaxies that 
the dust mass derived from extinction estimates is on aver- 
age ~ 29 x lower than that derived from FIR data. If this 
also holds for UGC 7321, then our estimated dust mass 
from our models appears to correlate roughly as expected 
with the estimate from the IRAS data,, especially given the 
uncertainties in the IRAS estimate (see, e.g., Lisenfeld & 
Ferrara 1998). 

Let us now return to the issue of the total interst ella r gas 
mass implied by our models. As described in Sect. 2.1 , our 



models are parameterized in terms of a total opacity and 
assume gas-to-dust ratio of 140 (Kim, Martin, & Hendry 
1994). As emphasized earlier, although the exact choice 



of the gas-to-dust ratio has negligible effect on the extinc- 
tion and reddening parameters derived from our models 
(since the opacity of the gas is negligible) , it does however 
impact the total gas+dust mass we infer: for p g / p d =140, 
we derive from Eq. 6 Mism ~ 1-2 x 10 8 for our f e B=4.0 
model. 

H I spectral line measurements by MGvD99 yield a total 
integrated Hi mass for UGC 7321 of Mhi=IA X 10 9 M G . 
Devereux & Young (1990) found that typically only the 
H I gas contained within roughly half of the optical radius 
of spiral galaxies has appreciable dust associated with it. 
From the pencil beam map of MGvD99, we estimate this 
corresponds to ~50% of the total H I content of UGC 7321, 
or -5.5xl0 8 M Q . In addition, Matthews & Gao (2000) 
recently detected 12 CO l-»-0 emission from the central 
regions of UGC 7321, and estimated a molecular hydro- 
gen content within the central 2.5 kpc of the galaxy of 
~ 2 x 10 7 M.q using a standard Galactic CO-to-H2 con- 
version factor (but see below). Thus after correction for 
helium, the total gas content in the central regions of 
UGC 7321 is ~6.2 times more than the total ISM mass 
inferred from our models. This discrepency is significant 
even after accounting for the range of accepted values for 
the Galactic gas-to-dust ratio [where p ga s/ Pdust ~ 120- 140 
based on depletion and extinction measurements (Whit- 
tet 1992)]. Not surprisingly, the assumption of a Galactic 
gas-to-dust ratio appears to be invalid for UGC 7321. 

There is now considerable evidence that the gas-to-dust 
ratio varies from galaxy to galaxy (e.g., Bouchet et al. 
1985; Hunter et al. 1989; Dwek 1998; Lisenfeld & Ferrara 
1998), with position within galaxies (e.g., Stanimirovic et 
al. 2000), and as a function of metallicity (e.g., Bouchet et 
al. 1985; Lisenfel & Ferrara 1998). For example, Koorn- 
neef (1982) found the gas-to-dust ratio for the Large Mag- 
ellanic Cloud (LMC) to be ~4 times the Galactic value, 
while Bouchet et al. (1985) and Stanimirovic et al. (2000) 
found a value for the SMC of ~8-ll times Galactic (as- 
suming a Galactic gas-to-dust ratio of ^140). 

Although no oxygen abundance measurements are avail- 
able for UGC 7321, typical metallicity values for LSB 
galaxies tend to be <; 1/3 solar (e.g., McGaugh 1994; 
Rdnnback & Bergvall 1995). That we require a gas-to- 
dust ratio of ~6.2 times the Galactic value to reconcile 
our model gas masses with observed values of UGC 7321 
thus seems appropriate for a galaxy whose metallicity is 
likely to lie somewhere between SMC and LMC values. 
The global gas-to-dust ratio inferred from our extinction 
estimates and the total observed Hi content of UGC 7321 
(after correction for He) would then be ^1600. 

3.6.2. Implications of a Clumpy, Multi-Phase ISM in 

UGC 7321 

We have shown that a dusty ISM with a significant 
fraction of its total mass (~56%) in a the form clumpy 
component is needed to explain the observed properties of 
UGC 7321. We now consider the implications of such a 
finding. 

For the type of simple, two-phase ISM model that we 
have used in the present work, it is often assumed that the 
clumped component of the medium should represent a pre- 
dominantly molecular gas medium, while the diffuse com- 
ponent represents diffuse atomic gas (e.g., Bianchi et al. 
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2000). In the case of UGC 7321, if we make such assump- 
tions, then we would predict a molecular gas to atomic 
gas ratio somewhat higher than the typical upper limits 
of <;40% found for the central regions of Sd-Sdm galaxies 
by Devereux & Young (1990). Moreover, this is consid- 
erably higher than empirical measures for UGC 7321 by 
Matthews & Gao (2000), which suggest that only ~4% of 
the total gas content in the central regions of UGC 7321 is 
molecular (assuming a Galactic CO-to-H 2 conversion fac- 
tor). 

There are several possible explanations for this trend. 
Firstly, the association of the atomic gas solely to a diffuse 
component, and the molecular gas solely to the clumped 
component in a gross oversimplification of the ISM of any 
galaxy (e.g., Turner 1988; Elmcgrccn 1993). This may be 
particularly true in the interstellar environments of LSB 
galaxies. The optical depths we infer for the "dark clouds" 
in UGC 7321 (A v £1.5; see also Gallagher et al. 2000) 
are consistent with the diffuse molecular clouds, and the 
lowest density Giant Molecular Clouds in the Milky Way 
(cf. Turner 1997). It has been suggested that such dark 
clouds in galaxies like the SMC may in fact contain sig- 
nificant atomic components (e.g., Dickey 1996; Kenncy 
1997). Even if such clouds are primarily molecular, their 
lower densities imply long timescales for the formation of 
CO molecules, and hence a reduction in the size of the 
CO-emitting regions within clouds of a given radius (e.g., 
Maloney 1990). In this case, CO measurements will sig- 
nificantly underestimate the molecular gas content of the 
clouds. Further uncertainties in the empirical H 2 estimate 
come from possible metallicity dependence on the CO-to- 
H2 conversion factor (e.g., Maloney 1990). A final possi- 
bility is that the gas-to-dust ratio differs between diffuse 
and clumpy ISM regions regions. This could occur, for 
example, if dust is accreted onto denser clouds, while dust 
in the diffuse medium is destroyed by shocks on timescales 
much shorter than new dust grains are created (e.g., Tie- 
lens 1998), although others have argued that gas and dust 
should be generally well mixed in the ISM (e.g., Cox & 
Mezger 1989). 

Regardless of the exact atomic versus molecular frac- 
tion in the clumpy ISM of UGC 7321, its very existence 
has interesting implications for furthering our understand- 
ing of the conditions in the ISM of LSB galaxies. Previ- 
ously, it has been suggested that it is likely that interstellar 
pressures in LSB galaxies would be too low to sustain a 
multi-phase structure in the diffuse, low metallicity envi- 
ronments in LSB galaxies (Mihos, Spaans, & McGaugh 
1999). Our data and models confirm that at least some 
LSB galaxies can maintain sufficient pressures to support 
a modest, multi-phase ISM. In addition, the existence of 
dust in ISM regions of enhanced density suggests that LSB 
galaxies have the envionments and the catalysts needed for 
molecular hydrogen formation. 

In general, a clumpy ISM should support higher H2 frac- 
tions than a nearly smooth ISM due to the presence of 
more efficient shielding of molecules from ionizing radia- 
tion fields (see also Spaans 1998). Thus even if the dark 
clouds in UGC 7321 contain significant atomic compo- 
nents, their inner regions almost certainly contain molecu- 
lar cores (as confirmed by the CO detection of Matthews & 
Gao) , hence eliminating the need to invoke star formation 



in LSB galaxies directly from atomic gas (cf. Schombert et 
al. 1990). This also raises the probability that other LSB 
galaxies could be detected in CO with sufficiently deep 
integrations (see also Matthews & Gao 2000). 

4. MODEL ROTATION CURVES 
4.1. Background 
4.1.1. The Importance of LSB Galaxy Rotation Curves 

Because of their relatively low visible matter contents 
relative to their inferred dynamical masses, LSB galaxies 
offer an important class of objects for unravelling the na- 
ture of dark matter in galaxies, as well as for testing dark 
matter alternatives, such as Modified Newtonian Dynam- 
ics (MOND; Milgrom 1983). An important step in this 
process is deriving accurate rotation curves. 

For some applications, systems that are nearly edge- 
on can be particularly desirable for rotation curve studies 
since inclination corrections to line-of-sight velocities and 
position angle uncertainties in these cases are small. H 11 
regions can also in general be observed to larger galacto- 
centric radii in edge-on galaxies than in less highly inclined 
systems, and edge-on galaxies are much easier to pick out 
and observe at higher redshifts. 

Unfortunately, the interpretation of observed rotation 
curves via mass-modelling techniques (e.g., Carignan 
1985) is strongly dependent on knowledge of the precise 
shape of the inner, rising portion of the rotation curve. 
This is often poorly determined observationally, particu- 
larly in edge-on galaxies. If rotation curves are derived 
from H 1 data, they often suffer from beam smearing (see, 
e.g., Swaters 1999 and references therein; van den Bosch 
et al. 2000) . High spatial resolution optical data can over- 
come this problem, but instead suffer from the vagaries 
and uncertainties of internal absorption in a given galaxy, 
and the debate over how accurately H II regions trace the 
true disk rotation (e.g., Davies 1991; Prada et al. 1994; 
McKeith et al. 1995; Prada 1996). Uncertainties due to 
projection effects and unknown gas distributions will also 
plague the interpretation of both optical and H I data in 
the edge-on case. 

Our Monte Carlo simulations provide a means to de- 
rive model optical rotation curves for edge-on galaxies and 
to test the severity of some of the influences mentioned 
above. Bosma et al. (1992) and Bosma (1995) have al- 
ready shown via rotation curve models and a comparison 
between H I-, CO-, and optically-derived rotation curves 
that even edge-on galaxies tend to be largely optically thin. 
Bosma et al. (1992), Byun (1993), and Bosma (1995) also 
showed that extremely high optical depths arc needed to 
produce an apparent linear rotation curve from an intrin- 
sically steeply rising one, even in edge-on and nearly edge- 
on galaxies. Here we reconfirm their results via Monte 
Carlo techniques, and also concentrate more specifically on 
exploring situations relevant to the interpretation of the 
rotation curve data for edge-on and nearly edge-on LSB 
galaxies, including UGC 7321. We also assess the promise 
of using such systems for detailed dynamical studies and 
mass modelling. 

4.2. Artificial Rotation Curves: the Method 

For the present study, we model galactic rotation curves 
by inputting an analytic galactic rotational law for our 
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model galaxy, and then adding together the Doppler shifts 
that arise between the emitted photons (presumed to arise 
from Ha-emitting "Hn regions") and the observer. We 
found that also including the Doppler shifts arising from 
the relative bulk motions of the scattering dust particles 
had negligible effect on the final rotation curves (see also 
Bosma et al. 1992). We assume that the dust opacity is 
constant over the narrow wavelength range that an ob- 
server would use to trace the galactic rotation. Because 
we are keeping track of wavelength in these simulations, 
when a photon exits the galaxy we project the photon into 
an image, and form "channel maps" by assigning an image 
to a particular wavelength range. 

The output of our radiation transfer code is a 3D data 
cube consisting of N two-dimensional channel maps. To 
form a rotation curve from our channel maps we take a 
strip across the major axis of each channel map and calcu- 
late the flux-weighted velocity at each position along the 
strip 

N N 

v(x)=Y,ViFi(x)/Y,Fi(x) , (9) 
i=i i=i 

where x is the position along the "slit" , is the velocity 
corresponding to the ith image channel, and Fi{x) is the 
flux at position x in the ith image channel. We note that 
this type of intensity- weighted determination of v(x) is 
technically correct only for non-edge-on galaxies (e.g., San- 
cisi & Allen 1979). However, it is still most commonly the 
method employed by optical observers of edge-on galax- 
ies (e.g., Goad & Roberts 1981; Karachentsev & Xu 1991; 
Makarov, Burenkov, & Tyurina 1999; Dalcanton & Bern- 
stein 2000), hence we adopt it here for ease of comparison 
of our models with data from the literature. The intrinsic 
spectral resolution of our models is 4.5 km s _1 pixel -1 . 
Before plotting our model rotation curves, we convolved 
them with Gaussians representing the expected turbulence 
of the ionized gas (assumed to have ct=10 km s _1 ) and a 
typical instrumental resolution (taken to be 30 km s _1 ). 
For simplicity, we begin by exploring smooth emissivity 
distri butio ns. The case of clumpy emissivity is explored in 
Sect. Oil. 



4.3. Artificial Rotation Curves: the Models 

As discussed above, the shape of the rotation curve ob- 
served for an edge-on galaxy at optical wavelengths may 
arise from a combination of effects, including: (1) the in- 
trinsic shape of the underlying rotation curve; (2) pro- 
jection effects; (3) dust extinction; (4) anomalously high 
extinction in the center only; (5) a complete absence of 
ionized gas in the central regions; (6) clumpy or irregu- 
larly distributed gas. Here explore the consequences of 
these effects for the observed rotation curves of edge-on 
galaxies similar to UGC 7321. 



4.3.1. 



and 



Effects of Internal Extinction, Projection, 
Intrinsic Rotation Curve Shape 

Let us first consider models that explore items (l)-(3) 
listed above, i.e., intrinsic rotation curve shape, projection 
effects, and wide-spread dust extinction. 

We take UGC 7321 as our template galaxy, adopting in 
all of our model rotation curves its inferred dust model 
distribution from Sect. 



2.4 



3.2 



and a maximum rota- 
tional velocity of 100 km s _1 (M"GvD99). For simplicity, 



we assume the scale height and scale length of the Ha 
corresponds to that of the stars. 

We begin with the simple case of a smooth dust emis- 
sivity distribution and consider edge-on optical depths of 
T"s,e=0,2.0, and 100. We compute our models for three 
inclinations: i = 85°, 88°, and 90°. As already shown by 
Bosma et al. (1992), the effects of items (l)-(3) above be- 
gin to become negligible at inclinations i <; 85°. 

For our input rotation laws, we choose four different 
rotation curve shapes. Three of the forms are based on 
the observed rotation curves of edge-on galaxies in the lit- 
erature having luminosities and peak rotational velocities 
similar to UGC 7321. The adopted models are (i) a very 
steeply rising curve, rising linearly to V max within 0.25 kpc 
and staying flat thereafter; (ii) a moderately steep curve, 
rising linearly to V rnax within the inner 1 kpc and flat- 
tening thereafter (e.g., UGC 3697, Goad & Roberts 1981); 
(Hi) an "arctan" rotation curve (see Courteau 1998) of the 
form v(R) = vq + -|u c arctan(i?) where R = (r — ro)/Vt, «o 
is the velocity center of rotation, ro is the spatial center of 
the galaxy, v c is an asymptotic velocity, and r t is a transi- 
tion radius between the rising and flat parts of the rotation 
curve (e.g., UGC 7321; Goad & Roberts 1981). For the 
arctan model we adopt «o=0, v c =100 km s _1 , ro=0, and 
r t = 1.3 kpc. Finally (iv) is a rotation curve that rises lin- 
early to 2.2h r and has only an essentially negligible flat 
portion within the stellar disk (see UGC 9242, Goad & 
Roberts 1981; FGC 1285, Dalcanton & Bernstein 2000). 

Our resulting models for these four different rotation 
curve shapes at three inclinations and three different op- 
tical depths are illustrated in Figure |l5|. In Figure 16 we 
also show Te^s—^-O models for input rotation curve forms 
(ii) and (Hi) with the Ha emission line data from Goad & 
Roberts (1981) overplotted. 

Several interesting trends are immediately evident in 
Figure |l5|. First, our model curves reconfirm the result of 
Bosma et al. (1992) that an observed rotation curve that 
is linear and rising to the outer regions of the stellar disk 
cannot be produced by any physically plausible amount 
of internal extinction. For i = 88°, we find that edge-on 
optical depths of tb.c ^100 are needed to produce such an 
observed rotation curve for any type of intrinsic rotation 
curve shape resembling models (i)-(iii)— i.e. for any shape 
other than an intrinsically slowly rising one. Moreover, 
it becomes clear from Figure [l5| that such a linear shape 
cannot be solely the result of deriving the rotational ve- 
locities via "intensity-weighted" procedures. We conclude 
that edge-on galaxies exhibiting linear, slowly rising rota- 
tion curves throughout their stellar disks (e.g., UGC 9242, 
Goad & Roberts 1981; FGC 1285, Dalcanton & Bernstein 
2000) must necessarily be galaxies with low central mat- 
ter densities and are likely to generally be optically thin 
systems. In the low optical depth regime appropriate for 
a typical LSB galaxy (fB, e ~4; see Sect. ??), the effects of 
dust extinction on the observed rotation curve are small 
enough to be safely ignored even in the case of i = 90°. 

At the opposite extreme, if a very steeply rising ro- 
tation curve that flattens appreciably in the inner few 
kpc is directly observed in a highly inclined galaxy (e.g., 
UGC 3697; Goad & Roberts 1981), Figure p] shows that 
one can infer that not only must the galaxy be relatively 
optically thin, but it must also be inclined at an angle 
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i <j 85° in order for its steep rise to be preserved in pro- 
jection. 

Figure [lq illustrates that either of the two intermedi- 
ate models[(n) or (Hi)] provides a roughly equally good 
match to the overall shape of the UGC 7321 rotation 
curve as measured by Goad & Roberts (1981), allowing 
for observational uncertainties and a slight lopsidedness in 
the true rotation curve. Thus we see that in these inter- 
mediate regimes, recovering the true, underlying rotation 
curve shape for a particular galaxy is less clear-cut. All 
of our models show that for the optical depths appropri- 
ate for a galaxy like UGC 7321, the effects of extinction 
can be safely ignored; instead, the greater ambiguity arises 
from the fact that rotation curves of quite different shapes, 
steepnesses and turnover radii can look very similar in pro- 
jection if the galaxy inclination angle i > 85°. Hence, this 
could create a much more serious problem for an observer 
attempting a dynamical model of a galaxy like UGC 7321. 
In the case where the inclination is not precisely known, 
the ambiguity becomes even more severe. 

4.3.2. Effects of Galaxy Centers Devoid of Gas 

It is sometimes suggested that the linear or slowly ris- 
ing rotation curves of some edge-on galaxies are due not 
to high optical depths throughout the disk, but rather 
a dearth of gaseous material in the inner regions of the 
galaxy, or the possibility that the H II regions (or H I gas) 
are confined primarily to a ring-like distribution. The for- 
mer case is essentially equivalent to a galaxy which for 
some reason may have a high optical depth near its cen- 
ter, but still may have an optically thin outer disk (cf. 
Davies 1991). Using our models, we explore the applicabil- 
ity of such models to UGC 7321, and to explaining the lin- 
ear, slowly rising rotation curves seen numerous late-type, 
edge-on galaxies (e.g., Goad & Roberts 1981; Karachent- 
sev & Xu 1991; Makarov et al. 1997; Makarov, Burenkov, 
& Tyurina 1999; Dalcanton & Bernstein 2000). 

We are immediately able to rule out a pure Hil region 
"ring" model for UGC 7321, and indeed, for most late- 
type, edge-on galaxies where Ha imaging data are avail- 
able. Such a ring morphology creates an Ha intensity 
distribution with peaks at the galactocentric radii corre- 
sponding to the radius of the projected ring; such a distri- 
bution is inconsistent with Ha imaging data available for 
UGC 7321 (MGvD99), for several other superthin galaxies 
(Ronnback & Bergvall 1995; Hoopes, Walterbos, & Rand 
1999; Matthews et al. unpublished) and for numerous less 
inclined LSB disk galaxies (e.g., Schombert et al. 1992; 
O'Neil et al. 1998). 

Although the ring model seems implausible, MGvD99 
did report a paucity of Ha emission near the center of the 
UGC 7321 disk, suggesting the possibility of a gas-poor 
or high optical depth region in the inner 0.5 kpc or so of 
the galaxy. For illustration purposes, here we test a some- 
what more extreme version of this model. In Figure 17, we 
show a model for i = 88° and rotation curve model (ii), 
but this time the model galaxy is completely devoid of Ha 
emission within a central hole of radius 1.5 kp c. All other 
parameters are as in the models in Figure [lq. 

It can be seen in Figure [l?] that the effect of this Ha- 
poor central disk region is to make the slope of the rotation 
curve slightly shallower than the case where the emission 
extends all the way to the center (cf. Figure O). However, 



this is still not sufficient to produce a linear rotation curve 
shape. The effect would of course increase somewhat if 
the hole were made larger, but such a model seems highly 
contrived; we are aware of no examples of Sd spiral galax- 
ies devoid of Ha emission over such a large fraction of the 
disk, nor should optical depth effects completely prevent 
Ha from being observable at small radii in normal late- 
type spirals, particularly LSBs. Hence this model should 
illustrate the upper limit of this effect for a realistic Sd 
galaxy. We conclude that the uncertainty of how far the 
emitting gas extends to the center of a galaxy will create 
additional small uncertainties in uncovering the the intrin- 
sic inner rotation curve shape for galaxies with i > 85°. 
However this effect cannot in itself produce an apparent 
linearly rising rotation curve for physically realistic galaxy 
models. 

4.3.3. The Effects of Clumpy of Sparsely Distributed Hu 

Regions 

Up until now, we have considered only smooth emis- 
sivity distributions for our rotation curve models. How- 
ever, the H II regions from which a rotation curve would 
be derived in a real galaxy are likely to have a clumpy 
and somewhat irregular distribution. Their distribution 
may be particularly sparse or irregular in the case of LSB 
galaxies (cf. Schombert et al. 1992; O'Neil et al. 1998). To 
explore how this may affect an observed rotation curve, we 
now compute a model where the Ha filling factor f fn a is 
only 0.5% and where C# Q =100. 

In Figure |l8| we show these "clumpy" models for the 
i = 88° case and input rotation curve form (ii). As ex- 
pected, the models now show a much more ragged appear- 
ance, with many of the fluctuations similar in amplitude 
to those in the real data. We see that an observed rota- 
tion curve with a significant number of "wiggles" is not 
necessarily a signature of strong patchy absorption, but 
could also occur in an extremely optically thin galaxy only 
sparsely populated with H II regions, as is not uncommon 
in the case of LSB galaxies. Our models suggest that these 
fluctuations may introduce some additional uncertainty in 
recovering the precise slope of the inner rotation curve of 
some LSB galaxies. 

5. SUMMARY 

We have presented some of the first 3D Monte Carlo 
simulations of the dusty ISM properties in the low optical 
depth regimes encountered in low surface brightness (LSB) 
spiral galaxies. We have demonstrated the power of using 
such techniques in combination with multiwavelength ob- 
servational data to constrain the amount and distribution 
of dust in such systems. Our realistic models fully take 
into account scattering and the effects of a clumpy, multi- 
phase ISM. 

Because late-type LSB galaxies are relatively optically 
thin, edge-on examples of these systems allow a unique op- 
portunity to explore the structure and distribution of their 
ISM and place important empirical constraints on dust 
models. We have shown that while dust contents in LSB 
galaxies are relatively low, these galaxies can contain mod- 
est amounts of dust and molecular material. Moreover, in 
at least some cases, a significant fraction of this dusty ma- 
terial is contained within a clumpy medium, confirming 
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that even LSB galaxies can maintain sufficient pressures 
to support a modest multi-phase ISM structure. 

From detailed models of the nearly edge-on LSB galaxy 
UGC 7321 (i = 88°), our Monte Carlo models indicate an 
edge-on optical depth f ei s ~4.0, and hence a total mass 
of dust is this galaxy of ~ 8.8 x 10 5 M Q . This is ~ 20x 
the warm dust content inferred from IRAS 100/xm mea- 
surements. Total 5-band extinction toward the center of 
this galaxy is estimated to be ~ 0.43 ± 0.2 mag. 

We infer that -50% of the dusty material in UGC 7321 
is contained in a clumpy medium. Based on the detec- 
tion of CO emission from UGC 7321 by Matthews & Gao 
(2000), it appears that some fraction of the clumped ma- 
terial is almost certainly molecular gas. However, the dark 
clouds in LSB galaxies may also contain significant frac- 
tions of atomic hydrogen gas. Together our models and 
direct Hi and CO measurements suggest that UGC 7321 
has a gas-to-dust ratio of at least 850 in the inner regions 
of its stellar disk, while the global gas-to-dust ratio for 
the entire galaxy (including H I gas extending beyond the 
stellar disk) is estimated to be —1600. 

In spite of the modest dust contents of LSB galaxies, we 
have shown that for the range of optical depths expected 
for such systems, dust extinction will have appreciable ef- 
fects on the observed total magnitudes and colors only 
when these galaxies are observed near edge-on. However, 
even then, the amount of reddening due to dust will be 
insufficient to explain the large radial color gradients now 
observed in a number of LSB galaxies (dc Blok et al. 1996; 
MGvD99; Matthews et al. 2000; Bell et al. 2000). Thus 
some LSB galaxies must have large intrinsic stellar popu- 
lation and/or metallicity gradients in their disks. 

We have shown that for realistic optical depths, dust 



has no appreciable effect on the rotation curves of edge-on 
spiral galaxies observed in the Ha emission line. Possi- 
ble holes in the central Ha distributions of such galaxies 
also cannot fully account for the observed slowly rising 
rotation curves. Thus the linearly, slowly rising rotation 
curves seen frequently in many late-type, edge-on or nearly 
edge-on spirals cannot be fully accounted for by dust and 
must be due to low central matter densities in these sys- 
tems. We have demonstrated that in general, projection 
effects will create a far greater uncertainty than optical 
depth effects in interpreting the precise intrinsic shapes of 
the rotation curves of edge-on galaxies. 
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Fig. 1. — F702W+F814W (R + I) composite image of a portion of the disk of UGC 7321, as imaged by the Wide Field and Planetary 
Camera 2 (WFPC2) aboard the Hubble Space Telescope. This is an 55" X 19" section of the disk, centered near r = 0. Resolution is ~0'.'l. 
White areas reveal the presence of clumpy, optically thick material in the disk of this LSB galaxy. 

Fig. 2. — Observed radial B — R color gradient along the major axis of UGC 7321. For ease of comparison with the model computations 
presented in this paper, the color has been shifted by a constant Ci such that (B — R) + c\ f«0 at the last reliably measured data points, near 
r =6 kpc. The data were extracted by averaging over a 12-pixel (2"3) wide strip, and smoothed by a factor of 15 for display purposes. For 
further details see MGvD99. 

Fig. 3. — Observed R — FF color gradient along the major axis of UGC 7321. The data were extracted by averaging along an 3" wide strip 
and then smoothed by a factor of 3 for display purposes. For ease of comparison with model calculations presented in this paper, the color 
has been shifted by a constant factor C2 such that (R, — FF) + C2 ~0 near the last reliable data points at r = ±3.5 kpc. 

Fig. 4. — R-band scattered light images generated from our Monte Carlo simulations of UGC 7321 at i = 88°. Each panel is 12.3 kpc 
is diameter, corresponding to the disk region over which our models were computed. The lefthand column shows models with smooth dust 
and emissivity distributions, and the righthand column shows mo de ls w ith a 2-phase clumpy ISM and stellar components. Further details on 
the input parameters for the models are summarized in Sect. p.l| & \'>.\ and Table 2. The four panels in each column represent four different 
edge-on B-band optical depths: r e s=8.0 (top); 4.0 (second row); 2.0 (third row); and 0.4 (bottom). 

Fig. 5. — Comparison of our best R- and H-band Monte Carlo model images of UGC 7321 with real data at these wavebands. Top: an 
i?-band image of UGC 7321 obtained with the WIYN telescope by MGvD99 and degraded in resolution to match the resolution of our models 
(~60 pc pixel -1 . Only the inner 12.3 kpc of the galaxy is shown. Row 2: an i?,-band model image with T e ^=4.0 and clumpy ISM and 
emissivity distributions. Third row: an H-band image of the inner 12.3 kpc of UGC 7321 obtained by MGvD99. Bottom: our model H-band 
image with r e g=A.O and clumpy ISM and emissivity distributions. 

Fig. 6. — Diff eren tial B — R. and R. — FF color profiles along the galaxy maj or axis for our family of i = 88° smooth ISM models (upper 
panels; see Sect. |3.l| ) and our clumpy ISM models (lower panels; see Sect. |3.3[ ). Shown are the color difference between models of 4 different 
edge-on optical depths r e g relative to a zero dust model, as a function of galactocentric distance: r e g = 8.0 (solid line); r e g = 4.0 (dotted 
line); r e g =2.0 (dashed line); r e g = 0.4 (dot-dash line). 

Fig. 7. — Illustration of the origin of the A(B — R) "saturation" effect for a slab model of uniform density and emissivity. Plotted are 
the color differences in B — R and R — FF between a constant density absorbing medium with no scattering and uniform emissivity, and a 
uniform dust-free emitting medium, as a function of B-band edge-on optical depth rg through the slab (see Equation 4). It can be seen that 
A(B — R) reaches a constant value at a considerably lower optical depth than A(R — FF). 

Fig. 8. — Differential FF — K color profiles along the galaxy major axis for our family of i = 88° clumpy models. Shown are the color 
difference between models of 4 different optical depths relative to a zero dust model, as a function of galactocentric distance. The 4 different 
line styles correspond to the same edge-on optical depths as in Figure <a. 

Fig. 9. — FF — K a color profile along the major axis of UGC 7321, derived from data from the Two Micron All-Sky Survey (2MASS). The 
data were extracted by averaging over a 3-pixel wide (~ 3") strip. The signal-to-noise of the data deteriorate significantly outside the central 
kpc of the galaxy (see Jarrett et al. 2000), but they are consistent with no significant FF — K s color gradient. For ease of comparison with 
model calculations presented in this paper, the observed color has been shifted by a constant factor C3 such that (FF — K s ) + C3 ~0. 

Fig. 10. — Dust-induced vertical B — R color gradients, extracted from our clumpy models parallel to the minor axis at three different 
disk positions. Shown are the color difference between models of 4 different edge-on optical depths T e g relative to a zero dust model, as a 
function of galactocentric distance: T Sl g = 8.0 (solid line); r e g = 4.0 (dotted line); T e g = 2.0 (dashed line); T e ,g = 0.4 (dot-dash line). The 
top panel shows the minor axis color profile; the middle panel shows the color at r=1.5 kpc, and the bottom panel along r = 3.0 kpc. 

Fig. 11. — Observed B — R vertical color profiles as a function of z (in kpc) for UGC 7321. The data were extracted over 15-pixel (~3"0) 
wide strips near the minor axis (top) and at r=3.0 kpc (bottom). For display purposes the data were smoothed by a factor of 3. For ease of 
comparison with the models in Figure llfjl the data have been shifted by a constant factor C4 such that (B — R) + C4 ssO near the last reliably 
measured data points, at r = ±0.4 kpc. 

Fig. 12. — B-, R-, iT-band normalized major axis intensity profiles for our clumpy models viewed at i = 88°, with T e s=0 (thick solid line); 
0.4 (dot-dash line); 2.0 (dashed line); 4.0 (dotted line); and 8.0 (solid line). 

Fig. 13. — Dust-induced B — R and R— FF color gradients as a function of galactocentric radius (in kpc) for smooth models models viewed 
at inclinations of i = 0°, 45°, & 90°. For ease of comparison, the i = 88° smooth models from Figure H are also reproduced here. Shown are 
the color difference between models of 4 different edge-on optical depths T e ,g relative to a zero dust model, as a function of galactocentric 
distance: r e t g = 8.0 (solid line); T e ,g = 4.0 (dotted line); r e g = 2.0 (dashed line); T B) g = 0.4 (dot-dash line). 

Fig. 14. — Dust-induced B — R and R — H color gradients as a fpiiction of galactocentric radius (in kpc) for clumpy models models viewed 
at inclinations of i = 0°, 45°, & 90°. Line styles are as in Figure Il3l For ease of comparison, the i = 88° clumpy models from Figure are 
also reproduced here. 

Fig. 15. — Rotation curves generated from our Monte Carlo models for 3 values of r e g (0, 2, & 100) and 3 viewing angles (i = 85 



& 90°) for the four analytic rotation curve models described in Sect. 1.3.1. Each column shows a different input model. The thick solid line 
in each panel shows the analytic form of the input rotation curve without any internal extinction. The thin solid line shows the "observed" 
intensity-weighted rotation curve for the r e g=0 (zero dust) model. The dotted line shows the predicted observed curve for a r e g=2 model, 
and the dashed line the predicted observed curve for a T e ^ = 100 model. In general, the r e g=0 and r e g=2 models are nearly indistinguishable. 

Fig. 16. — Rotation curves generated from our Monte Carlo models for r ei s=4.0 and a viewing angle of i = 88°. The thick solid line shows 
the form of the input rotation curve [model (ii) (top) and model (Hi) (bottom; see Text)], and the thin solid line shows the predicted observed 
rotation curve. Overplotted are datapoints (triangles) measured for UGC 7321 via Ha emission line spectroscopy by Goad & Roberts (1981). 

Fig. 17. — Model rotation curves generated for a galaxy seen at i = 88° and devoid of Ha emission over a central hole of radius 1.5 kpc. 
The thick solid line shows the form of the input rotation curve without any internal extinction. The thin solid line shows the "observed" 
intensity-weighted rotation curve for the r e ^g=0 (zero dust) model. The dotted line shows the predicted observed curve for a r e> g=2 model, 
and the dashed line the predicted observed curve for a T e g=100 model. Overplotted are data points (triangles) measured for UGC 7321 by 
Goad & Roberts (1981) via Ha emission line spectroscopy. 



16 



Fig. 18. — Model rotation curves generated for a galaxy seen at i = 88° and having a clumpy distribution of Hn regions, with a filling 
factor ff[ia=0.005. The thick solid line shows the analytic form of the input rotation curve without any internal extinction. The thin solid 
line shows the "observed" intensity-weighted rotation curve for the T Bt g=0 (zero dust) model. The dotted line shows the predicted observed 
curve for a T e ,.B=2 model, and the dashed line the predicted observed curve for a T ej s=100 model. Ovcrplotted are data points (triangles) 
measured for UGC 7321 via Ho emission line spectroscopy by Goad & Roberts (1981). 



Table 1. Properties of UGC 7321 



Hubble type: Sd IV 

V h : 407.2±2.2 (km s" 1 ) 
Distance: 10±| (Mpc) 

i : 88° 
a/b: 10.3 

D25. 5 : 5f6 (arcmin) 

^25.5: 16.3 (kpc) 

h r<R : 43"±4 w 2.1±0.2 (kpc) 

2'f9±0'f3 « 140 (pc) 
m s : 13.84±0.03 
B-i?: 0.99 
M B>i , : -17.05 
/xb(0): 21.6 (mag arcsec~ 2 ) 
Hba(0)'- 23.6 (mag arcsec~ 2 ) 
^_b: 27.6 (mag arcsec -2 ) 

M/ 20 : 233±7 (km s" 1 ) 
Mhi- 1.1x10 9 (A// q ) 
Mhi/Lb'- 1.1 (solar units) 

L F /a: 8x10 7 (L q ) 



Notes to Table 1. 

Adopted distance is from J. S. Gallagher (private 
comm.). The remaining parameters are taken from 
MGvD99 and Matthews (2000). Apparent magnitude 
tub, B — R color, and observed central surface bright- 
ness ^b(0) have been corrected for Galactic foreground 
extinction. Absolute B magnitude Mb, mean surface 
brightness Jib, and extrapolated central surface bright- 
ness ^B,i(0) were corrected for Galactic and internal 
extinction and projected to face-on values, as described 
in MGvD99. 
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Table 2. Model Input Parameters 



A 


k x (cm 2 /g) 


ax 


9\ 


C 


// 




//* 


B 


286 


0.66 


0.59 


20 


0.06 


2 


0.60 


R 


156 


0.53 


0.40 


20 


0.06 


2 


0.60 


H 


38 


0.42 


0.29 


20 


0.06 


2 


0.60 


K 


22 


0.36 


0.25 


20 


0.06 


2 


0.60 



Notes to Table 2. 

Explanation of columns: (1): waveband of the model cal- 
culation; (2) opacity of interstellar material; (3) albedo of 
interstellar grains; (4) phase function asymmetry param- 
eter; (5) density contrast between the clumpy and diffuse 
ISM (clumpy models only); (6) filling factor of the ISM 
clumps (clumpy models only); (7) density contrast between 
clumped and smooth emissivity (clumpy models only); (8) 
filling factor of the clumpy emission (clumpy models only) . 
For all models discussed in this paper: h r> * —2.1 kpc, 
h z ^—lA0 pc, h r>c t=15 kpc; h Zt d=70 pc, and both the dust 
and the stars were assumed to follow exponential vertical 
and radial distributions (see Sect. 2). 
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